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ABSTRACT 

The variability of active galactic nuclei (AGN) has long held the promise of shed- 
ding light on their detailed structure, and possibly other astrophysical phenonema. 
Different emission mechanisms lead to different patterns of variability in flux which 
are in principle easily distinguishable. Recent predictions for the expected spectrum of 
variations for various models are now in such a form that they can be compared with 
the observed statistical properties of AGN light curves from large scale monitoring 
programmes. In this paper, we use the results of a long term monitoring programme 
of a large sample of quasars and Seyfert galaxies, as well as individual light curves 
from the literature, to distinguish between the various model predictions. The results 
favour a model based on accretion disc instability over the starburst model where the 
variation comes from a succession of supernova bursts, but it also appears that much 
of the observed variation in quasars is due to gravitational microlensing. 
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1 INTRODUCTION 

Forty years after the discovery of quasars there is still much 
uncertainty about their structure. There is, arguably, broad 
agreement about the basic nature and overall arrangement 
of the various compone nts of a 'unified model' for active 
galactic nuclei (AGN) ( Antonucci 1991^ , but the details 
have proved very hard to tie down. In particular, the nature 



been carried out to measure the spectrum of variations in 
different wavebands. 

At present there are three basic models for explaining 
the observed AGN variability. The first involves instabilities 
in the accret ion disc, th e central engine powering the en- 
ergy output ( Rees 1984 ). The second postulates that AGN 
are powered by multiple supernova explosions or starbursts 



which result in stoch astic variations in brightness (Aretx 



are the subject ot mucn debate. Ir^art ot tne problem is tne 
paucity of observations which can effectively distinguish one 
model from another. Two types of constraint which are of 
particular interest are the spectral energy distribution and 
the observed variations in flux. Although one can deduce 
much about AGN st ructure from properties such as the op- 
tical/xray flux rat io Qcorgc fc Fabian 19911) or the so-called 
'big blue bump' (Gondhalekar et al., 1994|), such measure- 



of the [;entral engine and the radiative transter processes Terlevich 1994[). In the third approach, the observed 



ments are not sufficient to distinguish between competing 
models, let alone refine model parameters. 

Variations in fiux were detected in quasars shortly af- 
ter their discovery, and right from the start have played a 
pivotal role in constraining quasar morphology. The early 



detection of light fiuctuations on a timescale of months pro- 



vided a fundamental constraint on the size of the emitting 
region which has underlain all subsequent efforts to model 
the structure of quasars. In order to obtain a clearer picture 
of the emitting regions in AGN and the associated emis- 
sion mechanisms, a number of monitoring programmes have 



variations are not intrinsic to the AGN at all, but the re- 
sult of gravitational microlensing b y small compa ct bodies 
or MACHOs along the line of sight (Hawkins 199S). In fact, 
it seems quite likely that two or even all three of these pro- 
cesses are present at some level, and so the task of the ob- 
server is to disentangle them all in order to draw useful as- 
trophysical conclusions. There is also the possibility that dif- 
ferent mechanisms dominate in different luminosity regimes, 
and so in the analysis in this paper we shall divide AGN into 
two categories, quasars with Mb < —23 and Seyfert galaxies 
with A/fl > -23. 

Alth ough much effort has been put into monitoring 
quasars (Trcvcse et al., 1994; [Hook et al., 1994 ; Cristiani et 
al., 199^ ; p£^kmsT996 ), both individually and in samples 



of various sizes and different selection criteria, it seems fair 
to say that little light has so far been shed on the nature 
of quasars using these methods. There appears to be two 
main reasons for this. Firstly, despite the efforts which have 
been put into quasar monitoring programmes the data have 
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on the whole proved inadquate for measuring useful param- 
eters to characterise the variability. This is partly because 
samples of quasars have been too small, and partly because 
the run of data has been too short and inhomogeneous. As 
a consequence it has not proved possible to unambiguously 
define the fundamental properties of the flux variations of 
the quasars. 

The second reason that quasar monitoring programmes 
have not led to more progress in the understanding of AGN 
concerns the lack of firm predictions for variability from 
the various competing AGN models. A big step forward 
i n this area has recentl y been made with the publication 
(Kawaguchi et al., 1998) of a detailed model for AGN vari- 



ability from accretion disc instability. In this paper the au- 
thors make detailed statistical predictions for the spectrum 
of fluctuations for their model and also for the starburst 



model. These predictions are presented in a form that en- 



ables meaningful comparison with observations of AGN vari- 
ability. 



It is the purpose of this paper to use the best available 
observational data of the variations of AGN to distinguish 
between the various models of variability. To this end we 
describe a large scale monitoring programme of a sample 
containing some 600 quasars with regularly sampled light 
curves covering 24 years. We also use extensive monitoring 
data for Seyfert galaxies taken from the recent literature. 
Predictions of the models are published in the form of struc- 
ture functions and we analyse the observational data in the 
same way to provide quantitative comparisons. 



2 THE STRUCTURE FUNCTION 

In order to extract information from observations of AGN 
variability it is necessary to find ways of quantitively char- 
acterising the nature of the variations. There is a vast lit- 
erature on time series analysis, much of it concerned with 
extracting information from incomplete or inhomogeneous 
datasets. There are several functions which have been used 
for analysing AGN variations, including the structure func- 
tion, auto-correlation function and Fourier power spectrum. 
Although all of these functions contain very similar infor- 
mation, in practice they each have advantages in particular 
situations. For example, although for an infinite and com- 
plete run of data the auto-correlation function and power 
spectrum are essentially equivalent to each other, for finite 
runs of data there are significant differences. For long runs 
of evenly spaced data the power spectrum is to be preferred 
as it is on the whole easier to interpret and understand 
the errors. For short or inhomogeneous datasets the auto- 
correlation function provides a more stable measurement, 
but as the individual points are not independent of each 
other there can be difficulties with interpretation. The struc- 
ture function is very similar to the auto-correlation function 
and ha s been widely used i n the analysis of quasar light 
curves (iTrevese et al., 1994 |Hook et al., 1994| ; [Cristiani et| 
al., 1996 ; Arctxaga et a l., 1997) and microlensing statistics 



( Wyithc fc Turner 2001 ). The function of choice will depend 
upon a number of factors, not least of which is the form in 
which model predictions have been made. 

The pre sent paper has been la rgely prompted by the 
publication (Kawaguchi et al., 1998) of quantitative predic- 



tions for the statistics of AGN variability, and these are pre- 
sented in the form of structure functions. Accordingly, we 
shall proceed with the analysis of the observations in the 
same way. The structure function S may be defined by 



S{r) 



1 



Nir] 



El 



m{tj) - m(ti)]2 



where m{ti) is the magnitude measure at epoch ti, and the 
sum runs over the N{t) epochs for which tj — ti = t . The 
interpretation of the structure function in the sense of iden- 
tifying specific characteristics of the variation is not usu- 
ally feasible. However, particular models of variability can 
be shown to produce structure functions with measurable 
parametric forms. Although in some cases efforts h ave been 



made to predict the shape of the structure function ( Cid Fer- 



nandes et al., 1997), it is perhaps more useful to generate 



struc ture functions from simulated data (Kawaguchi et al 



199S), and this is the approach underlying the present pa- 
per. In addition to the standard structure function defined 
above, we shall for the purpose of measuring asymmetries 
also make use of two modified structure functions 5*+ and 
S-. These are defined as for 5* except that for S+ the inte- 
gration only includes pairs of magnitudes for which the flux 
becomes brighter, and for for which it becomes fainter. 



3 MODELS OF AGN VARIABILITY 

The mechanism behind AGN variability has been the sub- 
ject of much debate since their first discovery. At present 
there are three broad approaches to explaining the observed 
variations. The most favoured model has invoked instabili- 
ties in the accretion disc surrounding the central black hole, 
but a contending idea is that the variations are caused by 
some intermittent sequence of discrete events such as su- 
pernova bursts. A third possibility is that we are seeing the 
effects of gravitational microlensing by a population of small 
compact bodies along the line of sight. Each of these models 
has good arguments in its favour, but also fails to explain 
some aspects of AGN variability, and it seems likely that all 
three processes contribute at some level to the observed light 
curves with perhaps one mechanism dominating in any par- 
ticular regime. Until recently, there has been little attempt 
in the literature to make model predictions which can easily 
be tested against observations, but that is now changing. In 
the remainder of this section we review the current position 
for the three basic models mentioned above. 



3.1 Disk instability model 

The description of the accretion disc model as proposed by 
Rees (1984) contains much discussion about timescales of 
variability, but makes no predictions about the spectrum of 
variations. Since then there has been much work on variabil- 
ity in accretion discs (Wallinder et al., 1992), but the relation 



between the various types of instability and changes in emit- 
ted flux have proved hard to tie down and turn into specific 
testable predictions for an observed spectrum of variations. 
In fact to the extent that predictions are possible, they would 
appear to be at variance with much of the observed variabil- 
ity. 
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Figure 1. Simulated microlensing light curves from Fig. 2(a) of Lewis et al. (1993). 
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Figure 2. Structure functions for simulated microlensing light curves for model b of Schneider & Weiss (1987) (left hand panel), and 
from Fig. 2(b) of Lewis et al. (1993) (right hand panel). 



A major step forward has been made in a recent pa- 
per by Kawaguchi et al. (1998) who develop the cellular- 
automaton model for disc instability of Mineshige et al. 
(1994) to produce quantitative testable predictions for the 
shape of the structure function for the resulting photomet- 



ric variations of the disc. The basic idea of the cellular- 



automaton model is that as matter flows through an ac- 
cretion disc it causes instabilities which produce a spectrum 



of avalanches of different sizes, which in turn manifest them- 



selves as variations in emitted flux. The timescale of the light 
variation c orresponds to the accret ion timescale Tacc, which 
is defined (Kawaguchi et al., 1998) as 



160 



3/2 



M 
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where Tg is the Schwarzschild radius. For plausible values of 
the black hole mass M this implies a timescale of days or 
less for emission from the innermost stable orbit. However, 
observations of qua sar light curves imply timescales of at 



least several year s (Hook et al., 1994; Hawkins 199C; |Cris- 



tiani et al., 1996) which for the cellular-automaton model 



would mean emission from a characterist ic radial distance 
of~ lOOOr-g. In the notation of the paper (Kawaguchi et al., 
199^) this would imply Tin/ Tout ~ 1000. The authors use 
Monte Carlo simulations of the model for several parame- 
ter sets to produce simulated light curves from which they 
can evaluate structure functions. The structure functions 
have the form of a power law which flattens at the timescale 
Tacc- It may be seen from their Table 2 that the logarithmic 
slope of their structure function is not sensitive to the ratio 
rin/rout or to their other input parameter m' /m, the ratio of 
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Figure 3. Structure function for the simulated microlensing light 
curve from Fig. 2(b) of Lewis et al. (1993). 



diffusion mass to inflow mass, but lies close to 0.44 ± 0.03. 
This is a stable figure which can thus be compared with 
observations. 



3.2 Starburst model 

The idea that AGN variability might be caused by a series of 
discrete outbursts such as supernova explosions superposed 
as a Poisson process, has a long history. It was however 
recognised early on that there was a fundamental problem 
with the idea which had come to be known as the 'Christmas 



tree' model (Pica & Smith 1983). The model predicts a rela- 
tion between the luminosity and amplitude of quasars which 
is not observed, and so was ruled out early on as a serious 
contender for explaining their structure and variability. De- 
spite this the model has some attractive features, especially 
for explaining low luminosity AGN such as Seyfert galaxies 
where the nuclear luminosity could plausibly be explained by 
a few supernovae per year, which could at the same time ac- 
count for the observed amplitude of variability. For quasars 
of even moderate luminosity the required supernova rate of 
around one per day is such that the amplitude of variation 
would be far too small to be consistent with observations. 
Notwithstanding this, Terlevich et al. (1992) made a good 
case that the spectroscopic properties of AGN could be ex- 
plained as originating from supernovae with their 'starburst' 
model. The case that the variability of quasars can be ex- 
plained ( Aretxaga fc Terlevich 1994 ) is not convincing, but 
it cannot be dismissed as a possible dominating feature of 
Seyfert galaxies. 

One positive feature of the starburst model is the fea- 
sibility of predicting the spectrum of variations which will 



be observed. Kawaguchi et al. (1998), in addition to the ac- 
cretion disc model, also produced model light curves and 
structure functions for starbusts. They used the formal- 
ism of Aretxaga & Terlevich (1994) to construct a model 
and perform Monte Carlo simulations. Their Fig. 3 shows 
model light curves for a quasar of moderate luminosity which 
well illustrates the fundamental problem with the starburst 
model as applied to quasars. The amplitude of variation is 
about 0.2 magnitudes whi ch is far less t han would be ex- 
pected for such an object ( Hawkins 200C ). 

The structure functions of the model starburst light 
curves have a power law form with a flattening at around 
100 days, depending on model parameters. The power law 
section of the structure function has a logarithmic slope of 
about 0.83 ± 0.08, much larger than for the accretion disk, 
and a potential discriminant between the two models. 



3.3 Microlensing model 

The third broad approach to explaining observed AGN vari- 
ability invokes gravitational microlensing. Unlike the two 
models already described, no claim can be made that mi- 
crolensing explains all variations seen in AGN. The hypoth- 
esis is that for quasars the observed long term variations 
(on a timescale of several years) are dominated by the ef- 
fects of microlensing. Here, the main interest centres on the 
microlensing bodies, although there is the potential to learn 
about the strucure of the central region of quasars from the 
microlensing process. 

The idea ( [Hawkins 19"93| ; [Hawkins 199^ ) is that the dark 
matter is made up of a large population of planetary mass 
compact bodies or MACHOs which behave as cold dark mat- 
ter and randomly cross the line of sight to any particular 
quasar. The optical depth to microlensing is such that on 
average a quasar's light is lensed by several bodies at any 
one time, resulting in a distribution of caustics which causes 
a complicated amplification pattern on a timescale of a few 
years. Low redshift Seyfert galaxies are too nearby for any 
significant probability of microlensing, and so any variation 
seen in them must be intrinsic. Microlensing is well-known 
to occur in multiply lensed quasar systems where some vari- 
ations are only seen in one image, and on long timescales 
appear to dominate over the fiuctuatio ns seen in all im ages 
which must be intrinsic to the quasar (Pelt et al., 1995). 

Over the last 20 years or so a number of groups have 
published comput er simulations of the l i ght curves produced 



by rnicrolensing (Kayscr ct al., 1986 Schneider & Weisi 



1987] ; [Lewis et al., 1993| ). Fig. [i[ shows data from Fig. 2(a) of 



Lewis et al. (1993) of magnitude versus Einstein radius Re, 
sampled to enable comparison with observed light curves. 
Structure functions were calculated for these data, and also 
for the light curves for model b from Schneider & Weiss 
(1987). The two simulations are for point sources, but as- 
sume different surface mass density k« which i s essentially 
equivalent to the optical depth to microlensing ( Schneider & 



Weis s 1987). The two structu re functions are plotted in th e 
left ( pchncider fc Weiss 1987| ) and right ( [Lewis et al., 1993| ) 
panels of Fig. ^ The structure functions have logarithmic 
slopes of 0.28 and 0.23 respectively over the linear part, and 
there appears to be more power for greater optical depth. 

In microlensing simulations, the source size has a very 
noticeable effect on the appearance of the light curves, which 
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Figure 4. Light curves from the survey of Hawkins (1996), illustrating typical characteristics for quasars. 
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Figure 5. Structure function for the light curves of a sample of 
401 quasars from the survey of Hawkins (1996). 



have a smoother more rounded structure as it becomes 
larger. Fig. |^ shows the structure function for data from 
Fig. 2(b) of Lewis et al. (1993). These data come from the 
same distribution of lenses as for their Fig. 2(b), but with 
a source size of 0.2 Re- The logarithmic slope of the lin- 
ear part of the structure function in Fig. ^ is 0.31. This is 
slightly more than for a point source, as is to be expected 
from the loss of high frequency components attributable to 
a resolved source, although there is only a small difference 
in amplitude at long timescales. 



3.4 Comparison of model predictions 

It is interesting that the structure functions for all three 
models discussed above have a similar morphology. There is 
a power law or logarithmically linear rise from the shortest 
timescales to an eventual long timescale break or turnover. 
In each model the power or amplitude of the structure func- 
tion is somewhat dependent on the choice of model pa- 
rameters, but broadly speaking appears to be less for disc 
instability than for microlensing, and very much less for 
starburst. The timescales are difficult to compare, as it is 
only the starburst model which makes predictions for a spe- 
cific timescale, linked to the supernova cooling timescale 
which is well known, a characteristic value being 280 days 
(Kawaguchi et al., 1998). The accretion disc model is defined 
in terms of a time step which is left as a free parameter. Its 
value is related to the accretion timescale, which in turn de- 
pends on the size of the emitting region of the accretion disc. 
The timescale of microlensing models depends on the Ein- 
stein radius of the lenses and their mean transverse velocity. 
Although the second parameter can be estimated with rea- 



sonable confidence. Re is directly related to the mass of the 
lenses which in general is completely unknown. 

It is remarkable that the logarithmic slopes of the linear 
part of the structure functions show only small dispersions 
for different choices of parameter within each of the three 
models described above, but that the means are well sepa- 
rated. The slopes are 0.83 ±0.08, 0.44 ±0.03 and 0.25 ±0.03 
for the starburst, accretion disc and microlensing models 
respectively, which makes for a good opportunity for dis- 
tingishing between the models by comparison with observa- 
tions. 



4 OBSERVED STRUCTURE FUNCTIONS 
4.1 Quasars 

Over the past ten years or so there have been several ma- 
jor projects to monitor the light variations of quasars on a 
timescale o f ten years or more (Trevese et al., 1994 ; Hook 
et al., 1994 [Hawkins 199^; pristiani et al., 1996|).^e data 



for these monitoring programmes have been presented in 
the form of structure functions or auto-correlation functions, 
and used to look for correlations of redshift and luminosity 
with various parameters related to amplitude and timescale. 
It seems fair to say that the search for a correlation with the 
amplitude or power of the sructure function has been a lot 
more successful than questions involving timescale of varia- 
tion. 

The problem with the existing data is that the structure 
functions are not sufficiently well defined to compare with 
or distinguish between model predictions. There are a num- 
ber of reasons for this. Perhaps the most serious difficulty is 
the smallness of the timespan over which the quasars have 
been monitored, which means that features in the struc- 
ture functions such as breaks or slopes are not defined with 
sufficient accuracy. A related problem concerns lack of reg- 
ularity in the observations. Unevenly sampled light curves 
are notoriously difficult to analyse, resulting in various types 
of aliassing and artefacts which can be very misleading in 
their interpretation. A different problem concerns sample 
sizes. The difficulty of monitoring a large number of quasars 
individually on a regular basis is such that most monitoring 
projects have been based on wide field surveys, especially 
Schmidt telescopes. Even so most samples of quasars have 
been limited by available redshifts to around 300 members 
or less. A final limitation has concerned the analysis proce- 
dures used, which have been tailored to the unevenly spaced 
data that is hitherto all that has been available. 

In this paper we present data which go a long way to rec- 
tifying the problems mentioned above. The light curves are 
part of a monitoring programme which has been in progress 
since 1975, based on repeated measures of a single UK 1.2m 
Schmidt telescope field (ESO/SERC Field 287) centred on 
21h 28m -45°. The field has been observed several times ev- 
ery year in the Bj passband (Kodak Illa-J emulsion with 
a Schott GG395 filter) since 1977, as well as on shorter 
timescales of months, weeks and days, as well as in other 
passbands, notably R and U. In this paper we concentrate 
on the measures in Bj which form an unbroken sequence 
of 24 yearly measures from 1977 to 2000. At least one deep 
exposure was obtained every year during this period, and 
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Figure 7. Structure functions for the light curve of NGC 5548 from Peterson et al. (1999) and references therein for 50 day intervals 
(loft panel) and 10 day intervals (right panel). 
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Figure 8. Light curves from the survey of Hawkins (1996), illustrating typical characteristics for low redshift Seyfert galaxies. 



in most years the measures are based on four or more ex- 
posures. The plates were measured on the COSMOS or Su- 
perCOSMOS automated measuring machines at the Uni- 
versity of Edinburgh to produce yearly magnitudes for the 
~ 200, 000 objects in the central 19 square degrees of Field 
287. The photometric error on a measurement from a sin- 
gle plate was approximately 0.08 mag, and for most years 
when four plates were available this reduced to 0.04 mag. 
More basic details of the survey including a discusion of 



mea surement errors have already been published (Hawkini 



1996; f: awkins 2000), and here we concentrate on updating 
aspects of the survey relevant to the present paper. 

The quasars in the field have been detected using a 
number of techniques including ultra-violet excess, red drop- 
out, variability, objective prism searches and radio surveys. 
There are estimated to about 1500 quasars in the field to 
a magnitude limit of Bj = 22 of which 610 have been con- 
firmed with redshifts. Of these, several complete samples 
hav e been compile d according to various well-defined crite- 



ria (HajAfkms 200LJ 



The light curves show a variety of features, especially on 
a timescale of several years or more. Fig. ^ illustrates some 



examples. For most years the plotted magnitudes are the 
mean of four measures, and the error bars are based on the 
average photometric errors and the number of plates avail- 
able in any particular year, which in most cases was four. 
The light curves show no obvious assymetries in time, or 
easily definable morphological characteristics. Perhaps the 
one thing that can be said is that there appears to be more 
power on longer timescales, and that more time is needed to 
be sure that a characteristic timescale has been found. 

Fig. ^ shows the structure function for the quasar sam- 
ple. In order to make a clear distinction between Seyfert 
galaxies which are discussed below, a limit of z > 0.5 and 
AIb < —23 was set for the quasars, yielding a sample of 401 
objects. The errors were derived by splitting the sample into 
sub-units and measuring the dispersion in the associated 
structure functions. The structure function in Fig. ^ has a 
near power law form, with a logarithmic slope of 0.20±0.01. 
No allowance has been made for any time dilation eff'ect since 



this will not afi^ect the power law index ( Kawaguchi et al. 



1996 ), and it is easy to verify that in a narrow redshift range 
the index remains the same, albeit with a somewhat larger 
uncertainty due to the smaller sample size. 
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Figure 9. Structure function for the light curves of a sample of 
45 Seyfert galaxies from the survey of Hawkins (1996). 



There is no indication of a turn-over in the quasar struc- 
ture function at long timescales, which confirms the impres- 
sion from the hght curves in Fig. ^ that the quasars have 
not been monitored for long enough to reach a character- 
istic timescale. Equally there is no sign of the observations 
becoming dominated by noise at short timescales. This again 
confirms an impression from Fig. ^that the coherent varia- 
tions are much larger than the errors on the observations. 



4.2 Seyfert galaxies 

Seyfert galaxies are well known to vary in brightness on 
a timescale of months or less which makes them relatively 
easy to monitor within the constraints of modern policies 
for the allocation of t elescope time. Howev er, the project 
to monitor NGC 5548 



eterson et al., 199E| ) has surpassed 



all other efi'orts. Fig. |6 (top panel) shows the whole light 
curve sampled every 50 days for the purpose of measuring 
the structure function. The bottom panel shows the light 
curve from a particularly well observed period in 1993/4, 
sampled every 10 days. Error bars are not plotted as the 
errors are approximately the same size as the points of the 
plots. 

The length and frequency of observation of this light 
curve make it ideal for evaluating the structure function. 
The observations were made with spectrographs on a num- 
ber of telescopes as spectra were essential to the investiga- 
tors' main programme. The photometric observations were 
obtained by defining a window in the spectrum close to the V 
band and integrating within it. For most of the observations 
the errors are about 2% (0.02 magnitudes), but sometimes 
rise to twice this, about as large as the points in Fig. |^. 



The structure functions for the two light curves in Fig. ^ 
are shown in Fig. |^ Errors were estimated by re-sampling 
the data in various ways and measuring the resulting disper- 
sion in the structure functions. The left hand panel shows a 
rise to longer timescales with a gradually decreasing slope, 
becoming flat at a timescale of about a year. The logarith- 
mic slope does not become truly linear, and so to investigate 
the behaviour on short timescales we refer to the right hand 
panel for data sampled every 10 days. Here there is a well 
defined power law relation with logarithmic slope 0.38±0.01 
with an eventual flattening at lo nger timescales. 

The survey (Hawkins 1996) described in the previous 



section contains, in addition to the quasars, a number of 
Seyfert galaxies. Here again, in order to make a clear dis- 
tinction with the quasar sample we define our Seyfert galaxy 
sample with Mb > — 23 and z < 0.3 The light curves for 
this sample of 45 members have a very different character 
to those of the quasars. Fig. ^ shows some typical examples, 
with errors and other details as for Fig. ^ On the whole the 
light curves have small amplitudes with variations charac- 
terised by small fluctuations on short timescales, and smooth 
long term gradients of apparently modest amplitude. 

The structure function for this sample is shown in Fig. ^ 
with errors calculated as for Fig. ^. Again the data for 
Seyfert galaxies show a very different character from that 
for the quasars. The amplitude is much smaller as would be 
expected from the appearance of the light curves. There is 
a hnear section of the plot which has a slope of 0.36 ± 0.02, 
close to that of NGC 5548, and much larger than the figure 
of 0.20±0.01 for the quasars. The structure function flattens 
at short timescales. The reason for this is not immediately 
apparent, but may be related to measurement errors. 



5 TIME ASYMMETRIES IN STRUCTURE 
FUNCTIONS 

5.1 Model predictions 

In addition to the slope of the structure function, Kawaguchi 
et al. (1998) discuss a further discriminant between models 
of AGN variabiliy using the functions S- and S+ defined in 
section 2 above. These functions provide a measure of the 
underlying asymmetry of the emission process as manifested 
in the light curves. It is a fortuitous circumstance that the 
three models for variability discussed above make qualita- 
tively different predictions for the relationship between S- 
and S+. 

The disc instability model as put forward by Kawaguchi 
et al. (1998) predicts a spectrum of 'avalanches' in the ac- 
cretion disc which involve a gradual brightening of the disc 
followed by a sudden drop in flux. This pattern will mani- 
fest itself as a statistical asymmetry in the light curves which 
can be detected by comparing S- and S+. In this case S- 
should be larger than S+ towards shorter timescales. 

By contrast, the starburst model is dominated by the 
light of supernovae which are well known to rise rapidly in 
brightness at the onset of the explosion after which they 
gradually fade. Again, this asymmetry will be apparent in 
the AGN light curve as a statistical asymmetry, but this time 
s uch that S+ is larger th an S- towards shorter timescales 
( Kawaguchi et al., 1998 ). In both cases the difference be- 



© 2001 RAS, MNRAS 000, 



10 M. R. S. Hawkins 



Quasars 



NGC 5548 



CV3 



-1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — [— 



we 



_l I I I I I I I I I I I I I I I I I L_ 



1-0.5 0.5 1 

log t (years) 



1.5 



CO 

d 

I 



oo 



OO I 



CO 

oo 

I 



\ \ 


' ' 1 


1 1 1 


1 1 1 1 1 1 1 ^ 1 


1 


















A 










A 










A 




- 






A 


- 








_ AA 
























- 








- 








A » 






























T 










A 










A 








z 


















A 

, , 1 


1 1 1 


1 1 1 1 1 1 1 1 1 





1-2 -1.5 -1 -0.5 
log t (years) 



Figure 10. Time asymmetrical structure functions for the data in Fig. ^ (left panel) Fig. ^ (right panel). Functions for increasing (5+) 
and decreasing (S-) brightness are shown by upward and downward pointing triangles respectively. 



tween S- and 5*+ depends on the model parameters, but 
should be measurable with sufficiently good data. 

In the case of microlensing, which is an essentially 
symmetrical process in the context of quasar variability, it 
should not be possible to determine which way time is run- 
ning by any statistical tests on the fight curves. 5*- and S+ 
should be indistinguishable and lie on top of each other. 



5.2 Observations 



The idea of using time asymmetry in quasar light curves as 
a diagnostic for quasar vari ability has already received som e 
attention in the literature (Hawkins 1996; Aretxaga 1997). 
Tfie conclusion of both these studies is that statistically the 
light curves show no detectable departure from symmetry. 
The data used to produce the structure functions in Figs. ^ 
and ^ can also be used to calculate the functions S- and 
5*+. Fig. |l^ shows these two functions plotted together for 
comparison. Error bars have been omitted in the interests 
of clarity, but are as for the data in Figs. ^ and |^. 

The S- and S+ functions for quasars in the left hand 
panel lie very close to eacfi other, and are indistinguish- 
able within the errors. This implies that the observed vari- 
ation is symmetrical or time reversible. The data for the 
Seyfert galaxy NGC 5548 on the other hand show strong 
and significant differences between the two functions. On 
short timescales S- and S+ both have a logarithmic slope 
of about 0.4 but S- is systematically larger than S+. This 
implies that there is an asymmetry in the light curve in 
t he sense that the flux increases more slowly than it falls 
(Kawaguchi et al., 1998). On a time scale longer than three 



months the two curves separate, which is probably due to 
the limited length of the data run with 10 day sampling. 



6 CONFRONTATION OF MODELS WITH 
OBSERVATIONS 

We are now ready to test the three models described in sec- 
tion 3, which we shall do separately for quasars and Seyfert 
galaxies. The model predictions for the slope of the struc- 
ture function are 0.83 ± 0.08, 0.44 ± 0.03 and 0.25 ± 0.03 
for the starburst, accretion disc and microlensing models re- 
spectively. These figures are given in Table 1, along with the 
observed slopes for comparison. The slope measured for the 
quasar sample was 0.20 ± 0.01 which clearly rules out the 
starburst model. It also appears to be a long way from the 
accretion disc model but is consistent with the predictions 
for microlensing. It is however worth noting that for resolved 
sources the microlensing structure function gets steeper, pre- 
sumably as high frequency detail is blurred out, and there 
will certainly come a point where it is inconsistent with the 
data. 

We can further compare the models by testing for time 
asymmetries in the light curves. For short timescales the 
starburst, accretion disc and microlensing models predict 
S- < S+, S- > S+ and 5*- = S+ respectively. In fact as we 
have seen, for quasars S- and S+ are effectively coincident, 
again favouring the microlensing model. 

With Seyfert galaxies we find a very different picture to 
quasars. The slope of the structure function for NGC 5548 
and for the Seyfert galaxy sample can be taken as 0.38±0.01 
which again appears to rule out the starburst model, al- 



© 2001 RAS, MNRAS 000, |l|-0 



Variability in active galactic nuclei: confrontation of models with observations 11 



Table 1. Structure Function Slopes 



accretion disc models bec ause of the testable pre dictions 
which are mow available (Kawaguchi et al., 1998). There 



Model Predictions 


Model 


Slope 


Starburst 


0.83 ± 0.08 


Disc Instability 


0.44 ±0.03 


Microlensing 


0.25 ± 0.03 


Observations 


AGN Class 


Slope 


Seyfert Galaxies 


0.36 ± 0.02 


Quasars 


0.20 ±0.01 



though this time it lies close to the prediction of the accre- 
tion disc model. The microlensing model also appears to be 
excluded for any plausible value of the source size. Again we 
can further test the models by looking for time asymmetries, 
and we have found that for Seyfert galaxies S- > S+. This 
is in agreement with the prediction for the accretion disc 
model, and confirms the result for the structure function 
slope. 

To summarise, for Seyfert galaxies the observations 
favour the accretion disc model, and for quasars the mi- 
crolensing model. It appears that the starburst model is 
ruled out for all AGN. 



7 DISCUSSION 

Although within the terms of reference of this paper the re- 
sults in the previous section appear rather conclusive, there 
are a number of caveats which must be made clear. The pre- 
dictions for all three models were based on a very limited 
coverage of the relevant parameter spaces, and although in 
each case the dispersion in the structure function slope did 
not appear to be large, and very much less than the dif- 
ferences between the mean predictions of the models, it is 
entirely possible that a more thorough exploration of the 
parameter spaces could alter the picture significantly. No 
doubt in due course such data will become available, and 
the results given here can be tightened accordingly. 

In this paper we have concentrated on the slope of the 
structure function because for all three variability mecha- 
nisms considered it appears to be a robust quantity with It- 
tle dependence on the choice of model parameters. There are 
other quantities such as the timescale of variability or the 
amplitude of the turnover in the structure function which 
could in principle also be used. The problem here is the 
strong dependence on model parameters, or sensitivity to 
the window function of the data. We therefore content our- 
selves for the moment with the structure function slope. 
However, having once decided which variability mechanism 
is in operation, the process can be reversed to estimate the 
values of model parameters. 



The cellular-au tomaton model for disc instability (Mi 



neshige et al., 1994 ) has been used here as a prototype for 



is of course much controversy over the details of emission 
from an AGN accretion disc, and it may well be that other 
accretion disc models make very different predictions for the 
logarithmic slope of the structure function. When testable 
model data become available it seems possible that the con- 
clusions of this paper will have to be modified. It also goes 
without saying that the predictions of most models can be 
stretched to accommodate data by a careful choice of pa- 
rameters, and it may be that the conclusions of this paper 
turn out to be less conclusive than they appear to be at 
present. 

Notwithstanding these caveats, the conclusions reached 
in section 6 appear to be well supported by the data. Taking 
the results at face value, we find that photometric variation 
in Seyfert galaxies is best explained by instabilities in an 
accretion disc, while for quasars the dominant mechanism 
for variability is microlensing. At first sight this might seem 
paradoxical, but in fact the two explanations sit quite well 
together. We can account for them by a model in which 
the fluctuations in the light from the accretion disc be- 
come smaller with increasing luminosity, but the effects of 
microlensing become more pronounced at higher redshift, 
which for quasars typically means higher luminosity. 

This model can be tested directly in the rather rare 
situations where quasar variations can unambiguously be 
separated into those which are intrinsic to the quasar and 
those which are not. The best example for this purpose is 
the gravitational lens system 0957-1-561 which is a luminous 
quasar split into two components by the lensing effect of 
an intervening galaxy. The two images have been monitored 
intensely over the last 20 years or so, primarily for the pur- 
pose of determining the time delay between the two compo- 
i ients and hence me asuring the value of Hubble's constant 
(Kundic et al., 1997). The time delay was measured from the 
displacement of features which occured in both light curves 
and are clearly intrinsic to the quasar. These features are 
on the whole of short duration (~ 100 days) and small am- 
plitude (~ 0.1 magnitudes) as may be seen from Fig. 4 of 
Kundic et al. (1997). 

It has been known for a long time that as well as these 
intrinsic variations, the two components of 09574-561 also 
vary independently on a timescale of sev eral years with a n 
amplitude greater than 0.3 magnitudes (Pelt et al., 1998). 
This mode of variation is generally accepted as microlens- 
ing by compact bodies along the line of sight to the quasar 
images. The difference between the two modes of variation 
is well illustrated in Fig. 1 of Kundic et al. (1997) where 
the small repeating features contrast with the larger long 
term variation due to microlensing. These observations sup- 
port the idea that for quasars intrinsic variations are small 
(certainly smaller than for Seyfert galaxies) and that mi- 
crolensing can easily dominate the observed variations. 



8 CONCLUSIONS 

In this paper we have tested the predictions of three models 
for variability in AGN against observations. We have com- 
pared the predicted logarithmic slope of the structure func- 
tion for starburst, accretion disc instability and microlens- 
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ing models with the observed structure function for Seyfert 
galaxies and quasars. We find that for Seyfert galaxies the 
starburst and microlensing models are not compatible with 
the observed data, but the accretion disc model is in good 
agreement with it. For quasars the starburst and accretion 
disc instability models are ruled out while the microlensing 
model agrees well with the observations. 

We conclude that the variations in AGN are best ex- 
plained by a model in which at low luminosities the observed 
variations are caused by instabilities in the accretion disc. 
With increasing luminosity this type of variation becomes 
smaller in amplitude and the observed variation becomes 
dominated by the effects of microlensing. This increase is 
associated with the larger optical depth to microlensing for 
most quasars. 
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